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ABSTRACT 

Despite a century of remarkable progress in understanding stellar interiors, 
we know surprisingly little about the inner workings of stars spinning near their 
critical limit. New interferometric imaging of these so-called "rapid rotators" 
combined with breakthroughs in asteroseismology promise to lift this veil and 
probe the strongly latitude-dependent photospheric characteristics and even re- 
veal the internal angular momentum distribution of these luminous objects. Here, 
we report the first high precision photometry on the low-amplitude 6 Scuti vari- 
able star Rasalhague (a Oph, A5IV, 2.18 Mq, ^ ~ 0.88) based on 30 continuous 
days of monitoring using the MOST satellite. We have identified 57±1 distinct 
pulsation modes above a stochastic granulation spectrum with a cutoff of ~26 
cycles per day. Remarkably, we have also discovered that the fast rotation pe- 
riod of 14.5 hours modulates low-frequency modes (1-10 day periods) that we 
identify as a rich family of g-modes (|m| up to 7). The spacing of the g-modes 
is surprisingly linear considering Coriolis forces are expected to strongly dis- 
tort the mode spectrum, suggesting we are seeing prograde "equatorial Kelvin" 
waves (modes i = m). We emphasize the unique aspects of Rasalhague motivat- 
ing future detailed asteroseismic modeling - a source with a precisely measured 
parallax distance, photospheric oblateness, latitude temperature structure, and 
whose low-mass companion provides an astrometric orbit for precise mass deter- 
minations. 
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Introduction 



One of the most enduring questions of astronomy is "How do stars work?" Progress 
over the last century has led to a robust scientific framework that explains the physics 
of internal stellar structure and also how stars evolve in time. This framework strives to 
include stars of all types, e.g. low metallicity stars of the early universe, massive stars that 
will become supernovae, low-mass brown dwarfs as well as stars like our Sun. As a basic 
rule, our understanding is strongest and best proven for Sun-like stars and gets shakier and 
less robust for stars very much different from the Sun. This paper will highlight one critical 
aspect of stellar structure and evolution, but one that hardly affects our Sun - that of stellar 
rotation. 

Stars in general are probably born with significant angular momentum, but most of 
them are low-mass (like the Sun) and are "slowed" down by their magnetic winds in the 
early parts of their lives. On the other hand, most intermediate- and high- mass stars (>1.5 
solar masses) have weak magnetic fields and do not live very long, thus are often observed to 
be rotating very quickly. "Rapid rotati on" is expected to chang e the star's luminosity and 
photospheric temperature distribution (IMaeder fc Mevnetl 120001) and also strongly modify 
the observed surface abundances of various elements (jPinsonneaultl 119971 ) . For the most 
massive stars, rotation will partially determine when the star becomes supernovae, thus 
having a crucial impact on how heavy elements get dispersed into the interstellar medium. 
For all its importance, the effects of rapid rotation are only vaguely understood and cannot 
be confidently included in our stellar evolutionary codes. More observations are critically 



needed to support the recent renaissance in theoretical efforts 



and 3-dimensional hydrodynamical computer simulations fe.g.. iLee fc Saio 



20031 : 1 Jackson et al.ll2005l : iRoxburghl 12009 : lRieutordll2006l : iReese et al.ll2008h 



usmg ana. 



ytic calculations 



1997: Townsend 



Recent breakthroughs in making images of nearby rapidly-rotating stars using optical in- 
terferometry have coin cided with new numerical efforts to simulate the internal stellar struc- 



ture of these stars (e.g. . Ivan Belle et al.ll200ll : lDomiciano de Souza et al.ll2003l : iMonnier et al. 



*Based on data from the MOST satellite, a Canadian Space Agency mission operated by Dynacon, Inc., 
the University of Toronto Institute of Aerospace Studies, and the University of British Columbia, with 
assistance from the University of Vienna, Austria. 
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20071 : IZhao et al.ll2009l ). Many of the original assumptions that had been adopted are under 
renewed scrutiny and simple assumptions such as rigid-body rotation and early prescriptions 
for "gravity darkening" do not seem consistent with new data. 

Our group recognized the profound advances possible by co mbining the new i mages of 
nearby rapid rotators with the asteroseismic constraints (see also lCunha et al.ll2007l . for gen- 
eral discussion of astereoseismology and int erferometry) that c an be revealed by the precision 
photometer onboard the MOST satellite (jWalker et al.l 120031 ) . We have initially identified 
two rapidly rotating stars observable by MOST with extensive published interferometric 
datasets: Altair and a Oph (Rasalhague). Since there is no analytic theory that can predict 
the effects of rotation on the pulsational modes of a star spinning at nearly 90% of breakup, 
we plan to use the known stellar and geometric parameters of our sample to constrain the 
asteroseismic models with unprecedented power. Essentially, we know everything about the 
external appearance of these stars (size, oblateness temperature gradients on the surface, 
viewing angle, rotation periods) from interferometry and thus we can use asteroseismic sig- 
nal for deducing the internal structure without the normally disastrous degeneracies one 
confronts. 

Here we report on our first results of this project, namely the oscillation spectrum of 
a Oph using the MOST observations as part of the first year NASA Guest Observer program. 
a Oph is classified as an A5IV star (iGray et al.ll200ll ) with an estimated mass of 2.18 M© 
and is rotating at 88% of br eakup as judged by the extreme oblateness of the photosphere 
observed by interferometry (jZhao et al.ll2009l ). Interferometric imaging also revealed that 
our viewing angle is nearly exactly equator-on, a perspective that strongly suppresses the 
photometric amplitudes of / — |m| odd-parity modes, simplifying mode identifications. An 
early K star orbi ts the prirnary s tar with a period of 7.9 years (based on astrometric and 
speckle data, e.g. lGatewoodll2005l ). which will allow the masses of the two stars to be deter- 
mined eventually to high precision. Adaptive optics imaging has recently resolved (Hinkley 
et al., 2010, private communication) the two stars and finds the companion to be more than 
100 X fainter than the primary at wavelengths of the MOST instrument and so we have 
assumed all the pulsation signatures we see are from the primary. 

In this paper we present the full photometric dataset for Rasalhague along with an 
updated stellar model based on interferometry data. We find evidence for both p-mode and 
g-mode oscillations, along with a remarkable rotational modulation signal. It was beyond 
the scope of this paper to develop a mode analysis; however, we will pursue this in a future 
paper. 
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2. Properties of a. Oph 



The physical properties of a Oph we re recently determined using long-basehne inter- 
ferometry data f rom the CHARA Array (Iten Brummelaar et al.l |2005[ ) using the imaging 
combiner MIRC jMonnier et all l2004j . booeh . Izhao et al.l fcoooh reported that a Oph was 
being viewed within a few degrees of edge-on and with a rotation rate ~ 88% of breakup 
assuming a standard Roche model (point mas s gravity, sohd-body rotation) and a gravity 
darkening coefficient of 0.25 ( von Zeipel 1924al lbl). 



Since this paper was published, we have slightly improved the calibration and re- 
analyzed the existing interferometry data. Most importantly, we corrected a methodological 
error that resulted in inaccurate age and mass when comparing to non-rotating stellar evo- 
lution tracks. Figure [1] shows our new est imate of the location of a Oph o n the Hertzsprung- 
Russell diagram showing Y2 isochrones ( lYi et al.ll2001t iKim et al.l 120021 ). We have summa- 
rized the new model fits and the stellar parameters in Table [TJ Please see IZhao et al.l (120091 ) 
and Che et al (2010, in preparation) for a detailed description of the model parameters and 
for details on our fitting methodology. 

With the new analysis, we still find a Oph to be rotating 88% of breakup (rotation rate 
of 1.65 c/d) and being viewed within a few degrees of edge-on. However, the derived mass 
is now a bit higher - 2.18 Mq- a result of the improved isochrone analysis. We will adopt 
these parameters in discussions throughout the rest of the paper. Note that the error bars 
in the table rely on fixing the gravity darkening coefficient /3 to 0.25. Evidence is emerging 
(Xiao Che, private communication) that /3 might be substantially lower than this, which 
could affect the derived stellar parameters. 



Observations 



3.1. Basic Description 



The photometry presented here was obtained by the Microvariability and Oscillations of 
Stars (MOST) satellite. MOST houses a CCD photometer fed by a 15-cm Maksutov telescope 
through a custom broadband optical filter (350 - 750nm). The satellite's Sun-synchronous 
820km polar orbit (period = 101.4min = 14.20 cycles/day) enables uninterrupted observa- 
tions of stars in its Continuous Viewing Zone ( -18°< dec < -| -36°) for up to 8 weeks. A 



pre-launch summary of the mission is given bv IWalker et al.l (120031 ) and on-orbit science 
operations are described by iMatthews et al.l (120041 ). 



The MOST photometry of the a Ophiuchi binary system covers 30 days during 2009 
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May 27 to 2009 Jun 26, and was obtained in Fabry Imaging mode, where the telescope pupil, 
illuminated by the star, is projected onto the CCD by a Fabry microlens as an annulus 
covering about 1300 pi xels. The da t a wer e reduced using the background decorrelation 



technique described by iReegen et al.l (|2006[ ). The reduction includes filtering for cosmic 
particle impacts, correction of the stray light background modulated with satellite orbital 
parameters, filtering of long-term trends, and removal of obvious outliers. 

The final heliocentric corrected time series provided by the MOST instrument team 
contained 81018 measurements with an exposure time and sampling of 30 seconds (based 
on thirty 1-sec exposures stacked onboard the spacecraft). The datastream was remarkably 
consistent and continuous, being on-target with data for about 94% of the a Oph observation. 
Instrument problems did cause six data gaps in excess of 30 minutes, with the longest being 
5 hours. 



3.2. Data reduction 

The first step our team applied in our data analysis was to remove local 5-sigma out- 
liers from the data. Because a similar procedure had already been applied by the MOST 
instrument team, this step had only a minor effect - reducing the dataset to 80002 indepen- 
dent data points, a data volume reduction of only 1.3%. These flux points are presented in 
Figure [2] in their entirety. 

The next step was to r emove the ha r moni cs of the 14.2 c/d orbital period. Similar 



to the method described by I Walker et al.l (120051 ) . we used a moving window of 10 days to 



create a high signal-to-noise mean light curve folded at the orbital period. We then remove 
this component from the original light curve. The long averaging window ensures we only 
remove frequencies within a few (Nyquist) resolution elements of the harmonics. We have 
inspected data both with and without this correction and found the only significant effect 
was to remove the harmonics from our final power spectra. 

Lastly, we "detrended" the lightcurve in order to filter out slow variations on long time 
scales. The "slow trend" light curve was constructed by fitting a line segment to the flux 
measurements (after removal of orbital modulation) within a moving 1-day window. The slow 
trends were subtracted from the lightcurve resulting in a final data set on which subsequent 
Fourier analysis was performed. The slow trend light curve can be seen in Figure [H At 
our request, spacecraft telemetry was inspected to see if any camera or telescope diagnostic 
showed variations in sync with the slow variations seen here - none were found. 



We have graphically summarized the steps in our data reduction in Figure [31 Here we 
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have zoomed up on a single day in our time series and show exphcitly the original data, 
our estimate of the residual orbital modulation, and the slow trend. We also show the final 
corrected light curve with and without a 5-minute averaging window to demonstrate the 
overal quality of our data. 

After more than 7 years in orbit, the MOST Team has a very robust understanding of 
the performance of the instrument, especially the point-to-point precision of the photometry 
as a function of target magnitude and the noise levels at various frequency ranges in Fourier 
space. Many of the pulsating variable stars detected in the MOST Guide Star sample (in 
the magnitude ranges 6 < V < 10) had a "discovery amplitude" (i.e., amplitude of the 
largest oscillation signal) of only about 1-2 mmag, with many other oscillation peaks ten 
times smaller and a noise lev el in Fourier spac e 50 times lower (e.g., the richly multi-periodic 
delta Scuti star HD 209775 iMatthewd 120071 ). Some examples of stars which characterise 
the photometric precision of the MOST instrument are: HD 146490 (a guide star used as 
a calibration by iKallinger et al.l l2008bl ; see their Figure 1), with V = 7.2 and point-to- 
point scatter of 0.12 mmag; tau Bootis (V = 4.5), an exoplanet-hosting star which shows 
clear modulation at the planet's orbital period with a peak-to-peak amplitude of 2 mmag 
(jWalker et al.l l2008l ) for which its mean instrumental magnitude is repeatable to within 1 
mmag in observations spaced several years apart; and Procyon (V = 0.35), with point-to- 
point precision of 0.14 mmag. In these cases, the noise fioor in the Fourier spectrum is well 
defined and at a level of 0.01 mmag or less. 



Another MOST target, epsilon Oph ( iBarban et al.l 120071 : IKallinger et al.ll2008a[ ). a cool 
giant (V = 3.24), is stable to 0.35 mmag at frequencies below 0.1 c/d. At these low fre- 
quencies, non-periodic noise intrinsic to the star is usually associated with granulation. The 
MOST Team has never seen instrumental artifacts with the amplitudes and coherence of the 
variations evident at low frequencies in the MOST light curve of alpha Oph presented here. 
For the results in this paper, the low-frequency signals have no direct implications, so we do 
not address them in further detail. However, they are likely intrinsic to the star and hence 
do not raise any alarms of possible unrecognised photometric artifacts at other frequency 
ranges in the alpha Oph time series. 



4. Analysis 

4.1. Amplitude Spectrum 

We have carried out an exhaustive Fourier analysis of the MOST photometry. Although 
the MOST data are remarkably complete and evenly-spaced compared to ground-based data. 
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we nonetheless employed a numerical F o urier Transform followed by a "CLEAN" ing step 



based on the algorithm of [Roberts et al.l ( 119871 ) in order to remove the small sidelobes from 



the strongest peaks in the power spectrum that arise due to the gaps in the temporal sam- 
pling. 

Figure H] shows the amplitude spectrum of our photometry in units of milli-magnitude 
for the range to 64 c/d. The power for frequencies below ~1 c/d was suppressed by the 
detrending step described in the previous section. In fact, without detrending, nearly all 
the frequencies below ~0.75 c/d show significant power above background noise levels. We 
note that all frequencies up to 1440 c/d were inspected and frequencies above 64 c/d are not 
shown here because no statistically significant peaks were found. 



4.2. Non-white "noise" 

It was clear from the power spectrum that the background noise level is higher for 
frequencies below 30 c/d compared to higher frequencies. In order to properly estimate 
the statistical significance of a "peak" it is crucial to establish the level of the frequency- 
dependent noise background. We have used a median filter with a 2 c/d window in order 
to crudely estimate the background noise power for a given frequency, and this background 
level is included in Figure Si We will discuss the physical nature of this non-white noise in 
§5.1j for now we simply treat it as a noise source that can produce spurious spikes in the 
amplitude spectrum. 



4.3. Identification of Distinct Modes 



Armed with the frequency-dependent background noise spectrum estimated in the last 
subsection, we can now place confidence limits on a given amplitude peak being a "real" 
distinct pulsation mode rather than having arisen by chance due to background fluctuations. 
We have chosen a statistical measure such that we only expect one of our identified modes in 
the range between 0.5 and 64 c/d to be a false positive. Since there are ~1800 independent 
frequencies that can be probed in our data in this spectral range, we chose only peaks in 
the power spectrum above 3.5-cr confidence, a criteria that should produce roughly 1 false 
positive assuming pure gaussian random noise. The 57 modes that survived this analysis are 
listed in Table [2] and are plotted in Figure HI We also confirmed the statistical significance of 
our peaks through a bootstrap analysis and we note that we would have erroneously identified 
>500 distinct modes if we had ignored the frequency- dependence of the noise background. 
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The strongest modes are 18.668 c/d (half-amplitude 0.66 mmag), the closely spaced pair 
at 16.124/16.174 c/d (0.35/0.41 mmag), and at 11.72 c/d (0.41 mmag). These frequencies 
are typical for the dominant p-mode oscillations in 6 Scuti stars and be ar some resemblanc e 



to the dominant frequencies observed for another rapid rotator Altair ( iBuzasi et al.ll2005l ). 
which had dominant peaks at approximately 15.7, 20.8, and 26.0 c/d. The highest frequency 
mode we have detected with high confidence is at 48.35 c/d (0.036 mmag). 

Inspection of Figure S] reveals many additional peaks that appear real - more than can 
be attributed to random fluctuations of the background. However, we have chosen not to 
report these peaks in Table [2] in order to maintain a highly pure set of frequencies with 
minimal false positives. In other words, we understand that many of the other peaks are real 
pulsational frequencies but are also sure that some of them are contaminated by background 
fluctuations. Future mode analysis might allow some of the 2- or 3-0" peaks to become bona- 
fide pulsation frequencies, but we have been conservative in our mode identifications for this 
paper. 

In addition, inspection of the amplitude spectrum at lower frequencies easily reveals 
increased power at harmonics of 1.7 c/d, close to the expected rotational frequency (1.65 
c/d, see We will discuss the physical origin of these frequencies in the next section. 



4.4. Time Variability 

With such a long dataset of high quality, we investigated the stability of the detected 
modes. Figure S] contains also the power spectrum split into four 1-week chunks and is 
depicted as a 2-D background image in each panel. One can see the lower fr equency modes. 



<10 c /d, tend to be variable - strong in some weeks and not present in others. iDegroote et al. 



( 120091 ) found amplitude variability in modes in some higher-mass objects with liftimes from 
one to a few days, substantially shorter than those here which tend to be approximately one 
week. Most of the higher frequency modes are more stable, although not all. For instance 
the 16.1 c/d dominant mode shows amplitude varying by a large amount, consistent with 
two equal modes with a frequency difference just barely resolved by our 30 day run. The 
lifetime of modes is another clue to their physical origin and will be discussed in the next 
section. 
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5. Discussion 

Here, we wish to discuss three aspects of this dataset: the detection of granulation noise 
<30 c/d, the p-mode spectrum, and lastly the discovery of g-modes and their rotational 
modulation. 



5.1. Granulation 



Poretti et al.l fl2009f ) observed the 5 Scuti star HD 50844 using COROT, identifying 

>1000 statistically-significant pulsational modes using standard Fourier analysis. iKallinger fc Matthews 
( 120101 ) argued it was physically unreasonable to have so many distinct modes and interpreted 
the higher power for frequencies lower than ~30 c/d as due to granulation noise. Indeed, 
these authors presented a compelling correlation of the granulation cutoff frequency with 
stellar mass and radius. We subscribe to the interpretation of Kallinger & Matthews, choos- 
ing to treat the higher average Fourier power at the lower frequencies as a kind of non-white 
background noise (see previous section §4.21 

In order to better estimate the granulation cutoff frequency we have plotted the median- 
filtered power (2 c/d window) in Figure |5l We identified an obvious artifact from the orbital 
harmonics at 42.6 c/d, and 56.8 c/d and likely present at 14.2 c/d and 28.4 c/d, probably 
due to imperfect background subtraction of scattered light. By modeling the harmonics 
beyond 35 c/d, we have crudely estimated the contamination at lower frequencies based 
on the isolated features at 42.6 and 56.8 c/d. Figure |5] shows our ad hoc model for the 
contributions from the orbital harmonics and also includes our corrected granulation noise 
spectrum. Here, the granulation cutoff frequency can be estimated to be 26±2 c/d. 



As explained in IKallinger fc Matthews! (j2010[ ). the granulation cutoff frequency should 
be related to the inverse of the sound crossing time of one pressure scale height in the 



outer convective layer. This leads to a dependence of t'granuiation oc 



M 

i?2 ri 



It is beyond 



the scope of this paper to compare our result with their more elaborate power-law fitting 
method; however, we can compare our measured granulation cut off with that predicted by 
the empirical relation discovered by IKallinger fc Matthews! ( 120101 . see their figure 3). Using 
a mass of 2.18 M©, effective radius of 2.7R0and T^s = 83007^ for a Oph, we find that 
M R^"^ 2^-0.5 ^0.25 (in solar units), leading to an expected z/granuiation ~ 300;uHz= 26c/d, 
consistent with our observations. 
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5.2. Rotationally-modulated g-modes 



Following treatment of iBuzasi et al.l (|2005[ ). we can estimate the frequency of the fun- 



damental radial (p-) mode using the relation 



For Q = 0.033 days (lBregerlll979l ) and using the stellar parameters from Tabled], we 
find /fun ~10.1 c/d (compared to 15.6 c/d for the 1.79 Mq star Altair). Thus, we will begin 
by assuming p-modes all have frequencies above or near to this cutoff. A full analysis of the 
p-modes with tentative mode identifications will be the subject of a future paper. Here, we 
focus primarily on the unexpected discovery of rotationally-modulated g-modes in a Oph. 

We showed in the last section (see also Figure [2]) that there was a ~1 mmag slow (week- 
timescale) variation seen in the photometry. We also have identified 15 strong modes with 
frequencies smaller than the fundamental radial mode that cluster around harmonics of 1.7 
c/d, close to the rotation frequency. This distinctive structuring suggests we are witnessing 
the rotational modulation of modes whose frequency in the co-rotating frame is small in 
comparison to the rotation frequency. 

The tabulated frequencies were extracted directly from peaks in the CLEANed Fourier 
Transform and we assigned frequency errors based simply on the total length of the dataset. 
We realize that higher precision on the frequency localization of pure modes can sometimes 



be achieved through a multi-component least-squares fit (e.g. IWalker et al.l 120051 ). however 
this procedure has convergence problems for very closed spaced frequencies as encountered 
here. The errors on amplitudes were assigned based on the local background noise as es- 
timated in §4.2. We caution that the formal amplitude errors may be underestimates due 
to non-stationary statistics of the background or due to amplitude variations in the modes 
themselves. The full time series is available upon request to researchers interested in ex- 
ploring alternative analysis schemes, but our straightforward approach adopted here leads 
to conservative errors in the frequencies and is sufficient for the analysis presented here. 

Gravity (g-) modes are expected to have long periods compared to p-modes, although 
it has been difficult to estimate this theoretically for rapid rotators. Indeed, our amplitude 
spectrum is consistent with a complex and time- variable set of g-modes with co-rotating fre- 
quencies ^0.1 c/d, corresponding to approximately 10-day periods. Normally the properties 
of these modes are quite difficult to measure due to the long time frames needed to see the 
modes travel around the star. However, in our case, the ~1.65 c/d rotational rate of this 
rapid rotator leads to strong rotational modulation. Note also that modes of North-South 
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asymmetry are invisible to us while the symmetric I — \m\ even-parity modes are strongly 
modulated due to our near equator-on viewing angle. 

Specifically we see small clusters of modes around frequencies ~1.8, 3.5, 5.4, 7.0, 8.6, 
10.4, 12.0 c/d. While some of the higher frequencies here might be p-modes, it seems that 
the strong periodicity at about 1.7 c/d suggests that each cluster represents a new non-radial 
family of modes with increasing m values. Stellar rotation causes lightcurve variations with 
frequencies proportional to m/rot, since the m parameter represents the number of nodes 
around the equator. With this interpretation, we are seeing non-radial modes corresponding 
to m = 1, 2, 3, 4, 5 and possible m = 6, 7. We believe this is the first time such a rich set of 
g-modes has been detected and partially identified around a rapid rotato r, although perhap s 



a related phenomenon may be at work around the late Be star HD 50209 ( iDiago et al.ll2009l ). 



The simultaneous appearance of g-modes an d p-modes would make a Oph a "hybrid 



7 Dor-5 Set pulsator, a class recently explored by iGrigahcene et al.l (120101 ) using first Kepler 
data. Indeed, these authors present a frequency spectrum of KIC9775454 that bears some 
similarity to the data presented here for a Oph, although our frequency spectrum does not 
show the characteristic "gap" between 5 and 10 c/d. Further comparison with this work 
is not possible until more is known about the rotational properties of the Kepler sample, 
especially since many hybrid pulsators are known to be slowly-rotating Am stars. 

To further investigate the periodic spacing of the g-modes of a Oph, we searched for 
families of modes that follow a linear frequency relationship of the form 

/ = /o + m/i (2) 

for integer m. We found two families, each comprising more than three modes, that can 
be well represented by this formula. Family A contains 5 modes; assuming the lowest- 
frequency mode corresponds to m = 1, a least-squares fit gives /o = 0.073 ± 0.017 c/d and 
/i = 1.7095 ± 0.0002 c/d. Family B contains 4 modes, with a fit /o = 0.219 ± 0.017c/d and 
fi = 1.7424 ± 0.0002 c/d. These modes are noted in Table [2] and included in Figure along 
with the observed amplitude spectrum for / < 16 c/d. There are other solutions that fit 
three modes or fewer, but in this treatment we consider just Family A and B, since the larger 
number of modes mean the fit is more reliable and less likely due to chance superpositions. 

There are a number of possible interpretations for the linear frequency relationship 
exhibited by the two mode families, which we shall consider in turn. To first order in the 
rotation frequency /rot, the observed frequency of a mode with radial o rder n, harmonic 



degree i and azimuthal order m is given by the well-known iLedouxl ( 1l95ll ) formula 



/ = fn,e + mf T-otil - Cn,e)- 



(3) 
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Here, fn,e is the frequency the mode would have in the absence of rotation, and is independent 
of m. The parentheses combine together the effects of the Doppler shift (in transforming 
from co-rotating to inertial reference frame) and the Coriohs force — the latter represented 
by the Cn/ term, which again is ind ependent of m. For high-order g modes, Cn,e approaches 



+ (e.g.,|AertselalJ|2010|). 



If we make the identifications /o = /„/ and /i = /rot(l — C'n/), then the Ledoux formula 
is apparently able to explain the linear frequency relationship of the two families, assuming 
each comprises modes having the same n and i but differing m. However, there are strong 
arguments against this conclusion. For both families the mode frequencies in the co-rotating 
frame, 

fc = f- mf,ot (4) 

are small (on the order of 0.1 — 0.3 c/d) compared to the rotation frequency /rot = 1.65 c/d 
(Tabled]). The spin parameter u = 2|/rot//c| of the modes must therefore be much l arger than 



unity, implying that their dynamics are dominated by the Coriolis force (see, e.g. iTownsend 



20031 ) . This represents a fundamental inconsistency with the use of the first-order Ledoux 



formula Q, which is valid only for < 1. 

A possible way around this inconsistency is to relax the assumption that lowest-frequency 
mode in each family corresponds to azumithal order m = 1. If we instead assume m > 3 
for these modes, then by eqn. (jl])0 we can achieve u < 1. However, this 'fix' itself runs 
into difficulty when we note that for both mode families /i > /rot, implying that the Cori- 
olis coefficients C„.^ are negative. This is incompatible with the above-mentioned limit 
Cn,e ~ [^(^ + 1)]"^ > (indeed, we are not aware of any physically plausible pulsation model 
that predicts negative Coriolis coefficients). 

Accordingly, we are led to abandon the Ledoux formula as an explanation for the linear 
frequency relationship of the two mode families. In searching for alternative interpretations, 
we note that the small magnitude of /o suggests we may be observing 7 dispersion-free low- 
frequency modes. In rotating stars there are two classes of mode that are dispersion-free: 
equatorial Kelvin modes (which are prograde) and Rossby modes (which are retrograde). 
Rossby modes tend not to generate light variations be- cause they are incompressible, so let 
us put them aside for the moment (although not forget about them completely). 

Focusing therefore on the Kelvin modes, these are prograde sectoral {i = m) g-modes 
which are confined to an equatorial waveguide by the action of the Coriolis force (see 



^We note that negative /c values in this equation indicate modes whose phase propagation in the co- 
rotating frame is retrograde. 
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Townsendl l2003l ) . A defining characteristic of tliese modes is tliat tlieir frequencies are di- 



rectly proportional to m. At high radial orders, these frequencies in the co-rotating frame 
can be approximated by 

fc ~ TTT^; (5) 



here, following iMuUaru ( 1l989l ). we define 

h = 

with N is the Brunt- Vaisala frequency, 
given by 

f m 



N 



dr, 



V27r2 Jq r 

The corresponding observed frequencies are then 



27r2n 



rot 



(7) 



For a group of modes all having the same ra, this expression predicts a linear frequency 
relationship with zero intercept. If n varies within the group, however, then there will be 
some scatter about a linear relationship, with a relative amplitude of ~ An/n. 

To examine whether the two mode families could be Kelvin modes, we re-fit their fre- 
quencies assuming a zero intercept (again identifying the lowest-frequency mode as m = 1). 
This gives fjm = 0.07 ±0.03 c/d for Family A, and fjm = 0.10 ±0.03 c/d for Family B. To 
interpret these values, we assume a relationshi p I, ^ 10'^ W (M/Ma;)/(R/Ra,Y Hz, derived 
from the n = 3 polytropic model considered by iMuUaru (119891 ). Using the stellar parameters 
from Table [H this yields Jg ~ 34 c/d. Hence, the approximate radial orders of the modes are 
derived as n 24 (Family A) and n ~ 17 (Family B), with corresponding scatter (based on 
the quoted uncertainties in the fc/^ fits) of An ^ 10 and An ^ 5, respectively. 

These values are consistent with the typical radial o rders of unstable g modes observed 
in th e more-massive Slowly Pulsating B (SPB) stars (e.g.. lCameron et al.ll2008l : iPamyatnykh 



19991 ). This lends strong support to the identification of the two mode families as high-order 
[n ^ 17 and n ^ 24) g modes transformed by the Coriolis force into equatorial Kelvin modes. 

As a last comment, according to Figure HI most of the g-modes seem to be time- variable 
within our one-month time frame, in contrast to the p-modes whose amplitudes are more 
typically constant within the errors. The g-modes tend to change on a time scale similar to 
the intrinsic mode period (~10 days) or a bit longer. This could be due to superposition of 
many closely spaced modes in what appears to be indeed a complex spectra. Alternatively, 
many excitable modes may be transferring energy between themselves, leading to a dynamic 
and time-variable amplitude spectrum. The ability to see so many independent i and m 
modes by using strong rotational modulation is invaluable to deriving strong constraints on 
the mode properties up to high order (here up to m = 7) and promises to reveal interior 
stellar structure with additional analysis. 
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6. Conclusions 

This first long stare by MOST at a rapidly rotating A star has led to number of a new 
results. We detect for the first time a rich p-mode spectrum consistent with low-amplitude 
(5-Scuti pulsations, and measure a granulation spectrum below 26±2 c/d. In total, we have 
identified 57±1 distinct modes below 50 c/d including a complex set of low-frequency modes 
that we identify as rotationally- modulated g- modes with (co- rotating) frequencies ~0.1 c/d. 
A mode analysis revealed linear relationships between the spacings of g-modes up to m = 7, 
an unexpected result for a star rotating at ~90% of breakup. This periodicity can be 
explained as due to dispersion-free equatorial Kelvin waves (prograde / = m modes) although 
some inconsistencies in our analysis demand follow-up study. Lastly, the long time-base has 
allowed us to study the mode lifetime, finding that most p-modes are stable while g-modes 
appear to live only a few times their intrinsic (co-rotating) periods. 

Understanding the effect of rapid rotation on stellar interiors is crucial to developing re- 
liable models for massive star evolution in general, a Oph is emerging as a crucial prototype 
object for challenging our models and to spur observational and theoretical progress. Addi- 
tional work is planned using adaptive optics to determine the mass of the star by measuring 
the 8-year visual orbit more precisely, using visible and infrared interferometry to strictly 
constrain possible differential rotation and gravity darkening laws, and using asteroseismol- 
ogy to follow-up on the new discoveries outlined here. 

We thank C. Matzner for discussion and insights. TK is supported by the Canadian 
Space Agency and the Austrian Science Fund. AFJM is grateful for financial aid from 
NSERC (Canada) and FQRNT (Quebec). We acknowledge support from the NASA MOST 

guest observer program NNX09AH29G and NSF AST-0707927. This research has made use 
of the SIMBAD database, operated at CDS, Strasbourg, France, and NASA's Astrophysics 
Data System (ADS) Bibliographic Services." 
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Table 1. Best-fit and physical parameters of a Oph from CHARA-MIRC interferometry 



Model parameters 


Inclination (degs) 


87.5 ± 0.6 


Position Angle (degs) 


-53.5 ± 1.7 


Tpoi (K) 


9384 ± 154 


Rpoi (mas) 


0.757 ± 0.004 




0.880 ± 0.026 


P 


0.25 (fixed) 


Derived physical parameters 


Tequ (K) 


7569 ± 124 


Regu (R©) 


2.858 ± 0.015 


Rpoi (R©) 


2.388 ± 0.013 




8336 ± 39 


True Luminosity (L©) 


31.3 ± 0.96 


Apparent Tf^ff{K) 


8047 


Apparent Luminosity (Lq) 


25.6 


V Magnitude^ 


2.087 


H Magnitude'' 


1.709 


V sini (km/s) 


239 ± 12 


Rotation rate (cycles/day) 


1.65 ± 0.04 


Mass (M©)'^ 


2.18 ± 0.02 


Age (Gyrs)'^ 


0.60 ± 0.02 


xi of various data 


Total xi 


0.70 


CP xi 


1.08 


Vis2 xi 


0.86 


T3amp xi 


0.16 


Physical parameters from literature 


[Fe/H]'^ 


-0.16 


Distance (pc)" 


14.68 



^^V magnitude fro m literature; 2.086 ± 0.003 
jPerrvman et al.|[l997l 'l. 

H magnitude from literature: 1.72 ± 0.18 llCutri et al.l 
l2003h . 

"^Based on the stellar evolution model 

llDemaraue et ahlliool) . 

' ^Erspamer fc NorthI l l2003h 

teatewoodl ll2005h 
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Table 2. Distinct Modes detected for a Oph. 



Mode 


Center Frequency 


Half-Amplitude'* 


Mode 


# 


(cycles/day) 


(millimag) 


Family^ 


1 


1.768±0.017 


0.188±0.035 


A 


2 


1.835±0.017 


0.135±0.035 




3 


1.902±0.017 


0.152±0.033 




4 


3.428±0.017 


0.080±0.024 




5 


3.495±0.017 


0.160±0.025 


A 


6 


3.545±0.017 


0.108±0.024 




7 


3.695±0.017 


0.268±0.024 


B 


8 


5.439±0.017 


0.140±0.028 


B 


9 


6.923±0.017 


0.175±0.029 


A 


10 


7.024±0.017 


0.101±0.029 




11 


7.182±0.017 


0.207±0.029 


B 


12 


8.075±0.017 


0.183±0.033 




13 


8.508±0.017 


0.292±0.030 




14 


8.617±0.017 


0.115±0.028 


A 


15 


8.825±0.017 


0.224±0.028 




16 


10.227±0.017 


0.091±0.023 




17 


10.469±0.017 


0.104±0.025 




18 


10.619±0.017 


0.243±0.028 




19 


11.720±0.017 


0.405±0.033 




20 


12.028±0.017 


0.228±0.037 


A 


21 


12.412±0.017 


0.312±0.039 


B 


22 


13.096±0.017 


0.144±0.035 




23 


16.124±0.017 


0.349±0.039 




24 


16.174±0.017 


0.411±0.039 




25 


17.183±0.017 


0.115±0.027 




26 


18.209±0.017 


0.091±0.025 




27 


18.668±0.017 


0.655±0.026 




28 


18.818±0.017 


0.109±0.027 




29 


19.252±0.017 


0.111±0.026 




30 


19.936±0.017 


0.109±0.026 




31 


20.228±0.017 


0.093±0.028 




32 


20.286±0.017 


0.151±0.027 




33 


20.420±0.017 


0.120±0.027 




34 


20.512±0.017 


0.210±0.027 




35 


21.713±0.017 


0.282±0.032 




36 


22.155±0.017 


0.117±0.035 




37 


22.205±0.017 


0.299±0.034 




38 


22.480±0.017 


0.146±0.037 




39 


23.631±0.017 


0.223±0.027 




40 


23.807±0.017 


0.092±0.027 




41 


24.582±0.017 


0.136±0.027 




42 


25.166±0.017 


0.144±0.026 




43 


25.250±0.017 


0.114±0.027 




44 


25.416±0.017 


0.272±0.027 




45 


25.633±0.017 


0.111±0.027 
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Table 2 — Continued 



Mode 


Center Frequency 


Half-Amplitude^ 


Mode 


# 


(cycles / day) 


(millimag) 


Family'' 


46 


27.001±0.017 


0.112±0.030 




47 


29.120±0.017 


0.105±0.029 




48 


29.304±0.017 


0.172±0.030 




49 


30.980±0.017 


0.072±0.022 




50 


31.130±0.017 


0.089±0.022 




51 


32.949±0.017 


0.060±0.014 




52 


34.392±0.017 


0.049±0.015 




53 


35.718±0.017 


0.149±0.015 




54 


35.877±0.017 


0.113±0.015 




55 


39.597±0.017 


0.149±0.014 




56 


43.292±0.017 


0.157±0.022 




57 


48.347±0.017 


0.036±0.010 





"Amplitude error derived from level of frequency-dependent 
noise background (see §4.2). 

''See t|5.2l for descriptions of g-mode families A and B. 
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Fig. 1. — Based on re-analysis of CHARA-MIRC observations originally presented in 



we present a corrected HR diagram for a Oph compared to isochrones (|Yi et al.ll2001 



Zhao et al 



(1200911 ■ 



Kim et al. 



2002; 



Demargue et al.ll2004[ ). The asterisk marks the "apparent" position of a Oph, appearing red and faint due 
to the edge-on viewing angle. The solid line shows the range of possible apparent locations depending on 
viewing angle with the error bars marking the true luminosity and true effective temperature. Lastly, the 
square marks the non-rotating equivalent position for a Oph to allow comparison to the isochrones. 
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Alpha Oph Lightcurve from MOST 
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Fig. 2. — This figure shows the MOST photometry of a Oph (Rasalhague) in an early stage of data 
processing. Each tiny point is a measurement over a 30-second period and the sohd line shows the 1-day 
moving average. 
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Fig. 3. — This figure shows each step in the data reduction for the lightcurvc of a Oph (Rasalhague). The 
top panel shows the original data for an example 1-day period along with the estimated orbital modulation 
and slow trend components. The bottom panel shows the data after removing the orbital modulation and 
the slow trend - this final reduced light curve was used for Fourier analysis in this paper. 
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Fig. 4. — Here we present our Fourier Analysis of the lightcurve of a Oph. The soHd hue shows the 
amphtude spectrum based on the fuh 30-day data set, with the amphtude labeled on the right-hand axis. 
Our estimate of the median background noise level is shown with a dashed line and 3.5-(T peaks are shown 
with vertical lines. Here we also show the week- by- week temporal variability as a background image of the 
power spectrum, shown here in logscale. There were no statically-significant peaks discovered above 50 c/d. 
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The changing amplitude spectrum of Alpha Oph 
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Fig. 5. — This figure shows the median power spectrum for the light curve of a Oph using a 2 c/ci window 
(dot-dashed line). We note the orbital artifacts at harmonics of 14.2 c/d, and we have crudely estimated the 
orbital artifact (dashed line). We also plot the corrected granulation spectrum (thick solid line), identifying 
a clear frequency break at 26±2 c/d. 
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Fig. 6. — Here we repeat the amplitude spectruni of the detrended hght curve for frequencies below 16 c/d 
to highlight the g-modes. We also include the identifications for 2 families of modes which are discussed in 
the text C ^OI) . 



